KDM5-mediated transcriptional activation of ribosomal protein genes alters translation efficiency to regulate mitochondrial metabolism in neurons.

Genes encoding the KDM5 family of transcriptional regulators are disrupted in individuals with intellectual disability (ID). To understand the link between KDM5 and ID, we characterized five Drosophila strains harboring missense alleles analogous to those observed in patients. These alleles disrupted neuroanatomical development, cognition and other behaviors, and displayed a transcriptional signature characterized by the downregulation of many ribosomal protein genes. A similar transcriptional profile was observed in KDM5C knockout iPSC-induced human glutamatergic neurons, suggesting an evolutionarily conserved role for KDM5 proteins in regulating this class of gene. In Drosophila, reducing KDM5 changed neuronal ribosome composition, lowered the translation efficiency of mRNAs required for mitochondrial function, and altered mitochondrial metabolism. These data highlight the cellular consequences of altered KDM5-regulated transcriptional programs that could contribute to cognitive and behavioral phenotypes. Moreover, they suggest that KDM5 may be part of a broader network of proteins that influence cognition by regulating protein synthesis.

Molecular and network disruptions in neurodevelopment uncovered by single cell transcriptomics analysis of CHD8 heterozygous cerebral organoids.

About 100 genes have been associated with significantly increased risks of autism spectrum disorders (ASD) with an estimate of ~1000 genes that may be involved. The new challenge now is to investigate the molecular and cellular functions of these genes during neural and brain development, and then even more challenging, to link the altered molecular and cellular phenotypes to the ASD clinical manifestations. In this study, we use single cell RNA-seq analysis to study one of the top risk gene, CHD8, in cerebral organoids, which models early neural development. We identify 21 cell clusters in the organoid samples, representing non-neuronal cells, neural progenitors, and early differentiating neurons at the start of neural cell fate commitment. Comparisons of the cells with one copy of the CHD8 knockout and their isogenic controls uncover thousands of differentially expressed genes, which are enriched with function related to neural and brain development, with genes and pathways previously implicated in ASD, but surprisingly not for Schizophrenia and intellectual disability risk genes. The comparisons also find cell composition changes, indicating potential altered neural differential trajectories upon CHD8 reduction. Moreover, we find that cell-cell communications are affected in the CHD8 knockout organoids, including the interactions between neural and glial cells. Taken together, our results provide new data for understanding CHD8 functions in the early stages of neural lineage development and interaction.

Proteomics and phosphoproteomics profiling in glutamatergic neurons and microglia in an iPSC model of Jansen de Vries Syndrome.

Background: Jansen de Vries Syndrome (JdVS) is a rare neurodevelopmental disorder (NDD) caused by gain-of-function (GOF) truncating mutations in PPM1D exons 5 or 6. PPM1D is a serine/threonine phosphatase that plays an important role in the DNA damage response (DDR) by negatively regulating TP53 (P53). JdVS-associated mutations lead to the formation of a truncated PPM1D protein that retains catalytic activity and has a GOF effect because of reduced degradation. Somatic PPM1D exons 5 and 6 truncating mutations are well-established factors in a number of cancers, due to excessive dephosphorylation and reduced function of P53 and other substrates involved in DDR. Children with JdVS have a variety of neurodevelopmental, psychiatric, and physical problems. In addition, a small fraction has acute neuropsychiatric decompensation apparently triggered by infection or severe non-infectious environmental stress factors. Methods: To understand the molecular basis of JdVS, we developed an induced pluripotent stem cell (iPSC) model system. iPSCs heterozygous for the truncating variant (PPM1D+/tr), were made from a patient, and control lines engineered using CRISPR-Cas9 gene editing. Proteomics and phosphoprotemics analyses were carried out on iPSC-derived glutamatergic neurons and microglia from three control and three PPM1D+/tr iPSC lines. We also analyzed the effect of the TLR4 agonist, lipopolysaccharide, to understand how activation of the innate immune system in microglia could account for acute behavioral decompensation. Results: One of the major findings was the downregulation of POGZ in unstimulated microglia. Since loss-of-function variants in the POGZ gene are well-known causes of autism spectrum disorder, the decrease in PPM1D+/tr microglia suggests this plays a role in the neurodevelopmental aspects of JdVS. In addition, neurons, baseline, and LPS-stimulated microglia show marked alterations in the expression of several E3 ubiquitin ligases, most notably UBR4, and regulators of innate immunity, chromatin structure, ErbB signaling, and splicing. In addition, pathway analysis points to overlap with neurodegenerative disorders. Limitations: Owing to the cost and labor-intensive nature of iPSC research, the sample size was small. Conclusions: Our findings provide insight into the molecular basis of JdVS and can be extrapolated to understand neuropsychiatric decompensation that occurs in subgroups of patients with ASD and other NDDs.

Transcriptome analysis of neural progenitor cells derived from Lowe syndrome induced pluripotent stem cells: identification of candidate genes for the neurodevelopmental and eye manifestations.

BACKGROUND: Lowe syndrome (LS) is caused by loss-of-function mutations in the X-linked gene OCRL, which codes for an inositol polyphosphate 5-phosphatase that plays a key role in endosome recycling, clathrin-coated pit formation, and actin polymerization. It is characterized by congenital cataracts, intellectual and developmental disability, and renal proximal tubular dysfunction. Patients are also at high risk for developing glaucoma and seizures. We recently developed induced pluripotent stem cell (iPSC) lines from three patients with LS who have hypomorphic variants affecting the 3' end of the gene, and their neurotypical brothers to serve as controls. METHODS: In this study, we used RNA sequencing (RNA-seq) to obtain transcriptome profiles in LS and control neural progenitor cells (NPCs). RESULTS: In a comparison of the patient and control NPCs (n = 3), we found 16 differentially expressed genes (DEGs) at the multiple test adjusted p value (padj) < 0.1, with nine at padj < 0.05. Using nominal p value < 0.05, 319 DEGs were detected. The relatively small number of DEGs could be due to the fact that OCRL is not a transcription factor per se, although it could have secondary effects on gene expression through several different mechanisms. Although the number of DEGs passing multiple test correction was small, those that were found are quite consistent with some of the known molecular effects of OCRL protein, and the clinical manifestations of LS. Furthermore, using gene set enrichment analysis (GSEA), we found that genes increased expression in the patient NPCs showed enrichments of several gene ontology (GO) terms (false discovery rate < 0.25): telencephalon development, pallium development, NPC proliferation, and cortex development, which are consistent with a condition characterized by intellectual disabilities and psychiatric manifestations. In addition, a significant enrichment among the nominal DEGs for genes implicated in autism spectrum disorder (ASD) was found (e.g., AFF2, DNER, DPP6, DPP10, RELN, CACNA1C), as well as several that are strong candidate genes for the development of eye problems found in LS, including glaucoma. The most notable example is EFEMP1, a well-known candidate gene for glaucoma and other eye pathologies. CONCLUSION: Overall, the RNA-seq findings present several candidate genes that could help explain the underlying basis for the neurodevelopmental and eye problems seen in boys with LS.

Modeling the neuropsychiatric manifestations of Lowe syndrome using induced pluripotent stem cells: defective F-actin polymerization and WAVE-1 expression in neuronal cells.

Background: Lowe syndrome (LS) is a rare genetic disorder caused by loss of function mutations in the X-linked gene, OCRL, which codes for inositol polyphosphate 5-phosphatase. LS is characterized by the triad of congenital cataracts, neurodevelopmental impairment (primarily intellectual and developmental disabilities [IDD]), and renal proximal tubular dysfunction. Studies carried out over the years have shown that hypomorphic mutations in OCRL adversely affect endosome recycling and actin polymerization in kidney cells and patient-derived fibroblasts. The renal problem has been traced to an impaired recycling of megalin, a multi-ligand receptor that plays a key role in the reuptake of lipoproteins, amino acids, vitamin-binding proteins, and hormones. However, the neurodevelopmental aspects of the disorder have been difficult to study because the mouse knockout (KO) model does not display LS-related phenotypes. Fortunately, the discovery of induced pluripotent stem (iPS) cells has provided an opportunity to grow patient-specific neurons, which can be used to model neurodevelopmental disorders in vitro, as demonstrated in the many studies that have been published in the past few years in autism spectrum disorders (ASD), schizophrenia (SZ), bipolar disorder (BD), and IDD. Methods: We now report the first findings in neurons and neural progenitor cells (NPCs) generated from iPS cells derived from patients with LS and their typically developing male siblings, as well as an isogenic line in which the OCRL gene has been incapacitated by a null mutation generated using CRISPR-Cas9 gene editing. Results: We show that neuronal cells derived from patient-specific iPS cells containing hypomorphic variants are deficient in their capacity to produce F-filamentous actin (F-actin) fibers. Abnormalities were also found in the expression of WAVE-1, a component of the WAVE regulatory complex (WRC) that regulates actin polymerization. Curiously, neuronal cells carrying the engineered OCRL null mutation, in which OCRL protein is not expressed, did not show similar defects in F-actin and WAVE-1 expression. This is similar to the apparent lack of a phenotype in the mouse Ocrl KO model, and suggests that in the complete absence of OCRL protein, as opposed to producing a dysfunctional protein, as seen with the hypomorphic variants, there is partial compensation for the F-actin/WAVE-1 regulating function of OCRL. Conclusions: Alterations in F-actin polymerization and WRC have been found in a number of genetic subgroups of IDD and ASD. Thus, LS, a very rare genetic condition, is linked to a more expansive family of genes responsible for neurodevelopmental disorders that have shared pathogenic features.

Characteristics of allelic gene expression in human brain cells from single-cell RNA-seq data analysis.

BACKGROUND: Monoallelic expression of autosomal genes has been implicated in human psychiatric disorders. However, there is a paucity of allelic expression studies in human brain cells at the single cell and genome wide levels. RESULTS: In this report, we reanalyzed a previously published single-cell RNA-seq dataset from several postmortem human brains and observed pervasive monoallelic expression in individual cells, largely in a random manner. Examining single nucleotide variants with a predicted functional disruption, we found that the "damaged" alleles were overall expressed in fewer brain cells than their counterparts, and at a lower level in cells where their expression was detected. We also identified many brain cell type-specific monoallelically expressed genes. Interestingly, many of these cell type-specific monoallelically expressed genes were enriched for functions important for those brain cell types. In addition, function analysis showed that genes displaying monoallelic expression and correlated expression across neuronal cells from different individual brains were implicated in the regulation of synaptic function. CONCLUSIONS: Our findings suggest that monoallelic gene expression is prevalent in human brain cells, which may play a role in generating cellular identity and neuronal diversity and thus increasing the complexity and diversity of brain cell functions.

Transcriptome analysis of microglia in a mouse model of Rett syndrome: differential expression of genes associated with microglia/macrophage activation and cellular stress.

BACKGROUND: Rett syndrome (RTT) is a severe, neurodevelopmental disorder primarily affecting girls, characterized by progressive loss of cognitive, social, and motor skills after a relatively brief period of typical development. It is usually due to de novo loss of function mutations in the X-linked gene, MeCP2, which codes for the gene expression and chromatin regulator, methyl-CpG binding protein 2. Although the behavioral phenotype appears to be primarily due to neuronal Mecp2 deficiency in mice, other cell types, including astrocytes and oligodendrocytes, also appear to contribute to some aspects of the RTT phenotype. In addition, microglia may also play a role. However, the effect of Mecp2 deficiency in microglia on RTT pathogenesis is controversial. METHODS: In the current study, we applied whole transcriptome analysis using RNA-seq to gain insight into molecular pathways in microglia that might be dysregulated during the transition, in female mice heterozygous for an Mecp2-null allele (Mecp2+/-; Het), from the pre-phenotypic (5 weeks) to the phenotypic phases (24 weeks). RESULTS: We found a significant overlap in differentially expressed genes (DEGs) with genes involved in regulating the extracellular matrix, and those that are activated or inhibited when macrophages and microglia are stimulated towards the M1 and M2 activation states. However, the M1- and M2-associated genes were different in the 5- and 24-week samples. In addition, a substantial decrease in the expression of nine members of the heat shock protein (HSP) family was found in the 5-week samples, but not at 24 weeks. CONCLUSIONS: These findings suggest that microglia from pre-phenotypic and phenotypic female mice are activated in a manner different from controls and that pre-phenotypic female mice may have alterations in their capacity to response to heat stress and other stressors that function through the HSP pathway.

CRISPR/Cas9-mediated heterozygous knockout of the autism gene CHD8 and characterization of its transcriptional networks in cerebral organoids derived from iPS cells.

BACKGROUND: CHD8 (chromodomain helicase DNA-binding protein 8), which codes for a member of the CHD family of ATP-dependent chromatin-remodeling factors, is one of the most commonly mutated genes in autism spectrum disorders (ASD) identified in exome-sequencing studies. Loss of function mutations in the gene have also been found in schizophrenia (SZ) and intellectual disabilities and influence cancer cell proliferation. We previously reported an RNA-seq analysis carried out on neural progenitor cells (NPCs) and monolayer neurons derived from induced pluripotent stem (iPS) cells that were heterozygous for CHD8 knockout (KO) alleles generated using CRISPR-Cas9 gene editing. A significant number of ASD and SZ candidate genes were among those that were differentially expressed in a comparison of heterozygous KO lines (CHD8+/-) vs isogenic controls (CHD8+/-), including the SZ and bipolar disorder (BD) candidate gene TCF4, which was markedly upregulated in CHD8+/- neuronal cells. METHODS: In the current study, RNA-seq was carried out on CHD8+/- and isogenic control (CHD8+/+) cerebral organoids, which are 3-dimensional structures derived from iPS cells that model the developing human telencephalon. RESULTS: TCF4 expression was, again, significantly upregulated. Pathway analysis carried out on differentially expressed genes (DEGs) revealed an enrichment of genes involved in neurogenesis, neuronal differentiation, forebrain development, Wnt/ß-catenin signaling, and axonal guidance, similar to our previous study on NPCs and monolayer neurons. There was also significant overlap in our CHD8+/- DEGs with those found in a transcriptome analysis carried out by another group using cerebral organoids derived from a family with idiopathic ASD. Remarkably, the top DEG in our respective studies was the non-coding RNA DLX6-AS1, which was markedly upregulated in both studies; DLX6-AS1 regulates the expression of members of the DLX (distal-less homeobox) gene family. DLX1 was also upregulated in both studies. DLX genes code for transcription factors that play a key role in GABAergic interneuron differentiation. Significant overlap was also found in a transcriptome study carried out by another group using iPS cell-derived neurons from patients with BD, a condition characterized by dysregulated WNT/ß-catenin signaling in a subgroup of affected individuals. CONCLUSIONS: Overall, the findings show that distinct ASD, SZ, and BD candidate genes converge on common molecular targets-an important consideration for developing novel therapeutics in genetically heterogeneous complex traits.

Integrative transcriptome network analysis of iPSC-derived neurons from schizophrenia and schizoaffective disorder patients with 22q11.2 deletion.

BACKGROUND: Individuals with 22q11.2 Deletion Syndrome (22q11.2 DS) are a specific high-risk group for developing schizophrenia (SZ), schizoaffective disorder (SAD) and autism spectrum disorders (ASD). Several genes in the deleted region have been implicated in the development of SZ, e.g., PRODH and DGCR8. However, the mechanistic connection between these genes and the neuropsychiatric phenotype remains unclear. To elucidate the molecular consequences of 22q11.2 deletion in early neural development, we carried out RNA-seq analysis to investigate gene expression in early differentiating human neurons derived from induced pluripotent stem cells (iPSCs) of 22q11.2 DS SZ and SAD patients. METHODS: Eight cases (ten iPSC-neuron samples in total including duplicate clones) and seven controls (nine in total including duplicate clones) were subjected to RNA sequencing. Using a systems level analysis, differentially expressed genes/gene-modules and pathway of interests were identified. Lastly, we related our findings from in vitro neuronal cultures to brain development by mapping differentially expressed genes to BrainSpan transcriptomes. RESULTS: We observed ~2-fold reduction in expression of almost all genes in the 22q11.2 region in SZ (37 genes reached p-value < 0.05, 36 of which reached a false discovery rate < 0.05). Outside of the deleted region, 745 genes showed significant differences in expression between SZ and control neurons (p < 0.05). Function enrichment and network analysis of the differentially expressed genes uncovered converging evidence on abnormal expression in key functional pathways, such as apoptosis, cell cycle and survival, and MAPK signaling in the SZ and SAD samples. By leveraging transcriptome profiles of normal human brain tissues across human development into adulthood, we showed that the differentially expressed genes converge on a sub-network mediated by CDC45 and the cell cycle, which would be disrupted by the 22q11.2 deletion during embryonic brain development, and another sub-network modulated by PRODH, which could contribute to disruption of brain function during adolescence. CONCLUSIONS: This study has provided evidence for disruption of potential molecular events in SZ patient with 22q11.2 deletion and related our findings from in vitro neuronal cultures to functional perturbations that can occur during brain development in SZ.

Reduced CYFIP1 in Human Neural Progenitors Results in Dysregulation of Schizophrenia and Epilepsy Gene Networks.

Deletions encompassing the BP1-2 region at 15q11.2 increase schizophrenia and epilepsy risk, but only some carriers have either disorder. To investigate the role of CYFIP1, a gene within the region, we performed knockdown experiments in human neural progenitors derived from donors with 2 copies of each gene at the BP1-2 locus. RNA-seq and cellular assays determined that knockdown of CYFIP1 compromised cytoskeletal remodeling. FMRP targets and postsynaptic density genes, each implicated in schizophrenia, were significantly overrepresented among differentially expressed genes (DEGs). Schizophrenia and/or epilepsy genes, but not those associated with randomly selected disorders, were likewise significantly overrepresented. Mirroring the variable expressivity seen in deletion carriers, marked between-line differences were observed for dysregulation of disease genes. Finally, a subset of DEGs showed a striking similarity to known epilepsy genes and represents novel disease candidates. Results support a role for CYFIP1 in disease and demonstrate that disease-related biological signatures are apparent prior to neuronal differentiation.

CRISPR/Cas9-mediated heterozygous knockout of the autism gene CHD8 and characterization of its transcriptional networks in neurodevelopment.

BACKGROUND: Disruptive mutation in the CHD8 gene is one of the top genetic risk factors in autism spectrum disorders (ASDs). Previous analyses of genome-wide CHD8 occupancy and reduced expression of CHD8 by shRNA knockdown in committed neural cells showed that CHD8 regulates multiple cell processes critical for neural functions, and its targets are enriched with ASD-associated genes. METHODS: To further understand the molecular links between CHD8 functions and ASD, we have applied the CRISPR/Cas9 technology to knockout one copy of CHD8 in induced pluripotent stem cells (iPSCs) to better mimic the loss-of-function status that would exist in the developing human embryo prior to neuronal differentiation. We then carried out transcriptomic and bioinformatic analyses of neural progenitors and neurons derived from the CHD8 mutant iPSCs. RESULTS: Transcriptome profiling revealed that CHD8 hemizygosity (CHD8 (+/-)) affected the expression of several thousands of genes in neural progenitors and early differentiating neurons. The differentially expressed genes were enriched for functions of neural development, ß-catenin/Wnt signaling, extracellular matrix, and skeletal system development. They also exhibited significant overlap with genes previously associated with autism and schizophrenia, as well as the downstream transcriptional targets of multiple genes implicated in autism. Providing important insight into how CHD8 mutations might give rise to macrocephaly, we found that seven of the twelve genes associated with human brain volume or head size by genome-wide association studies (e.g., HGMA2) were dysregulated in CHD8 (+/-) neural progenitors or neurons. CONCLUSIONS: We have established a renewable source of CHD8 (+/-) iPSC lines that would be valuable for investigating the molecular and cellular functions of CHD8. Transcriptomic profiling showed that CHD8 regulates multiple genes implicated in ASD pathogenesis and genes associated with brain volume.

MicroRNA Profiling of Neurons Generated Using Induced Pluripotent Stem Cells Derived from Patients with Schizophrenia and Schizoaffective Disorder, and 22q11.2 Del.

We are using induced pluripotent stem cell (iPSC) technology to study neuropsychiatric disorders associated with 22q11.2 microdeletions (del), the most common known schizophrenia (SZ)-associated genetic factor. Several genes in the region have been implicated; a promising candidate is DGCR8, which codes for a protein involved in microRNA (miRNA) biogenesis. We carried out miRNA expression profiling (miRNA-seq) on neurons generated from iPSCs derived from controls and SZ patients with 22q11.2 del. Using thresholds of p<0.01 for nominal significance and 1.5-fold differences in expression, 45 differentially expressed miRNAs were detected (13 lower in SZ and 32 higher). Of these, 6 were significantly down-regulated in patients after correcting for genome wide significance (FDR<0.05), including 4 miRNAs that map to the 22q11.2 del region. In addition, a nominally significant increase in the expression of several miRNAs was found in the 22q11.2 neurons that were previously found to be differentially expressed in autopsy samples and peripheral blood in SZ and autism spectrum disorders (e.g., miR-34, miR-4449, miR-146b-3p, and miR-23a-5p). Pathway and function analysis of predicted mRNA targets of the differentially expressed miRNAs showed enrichment for genes involved in neurological disease and psychological disorders for both up and down regulated miRNAs. Our findings suggest that: i. neurons with 22q11.2 del recapitulate the miRNA expression patterns expected of 22q11.2 haploinsufficiency, ii. differentially expressed miRNAs previously identified using autopsy samples and peripheral cells, both of which have significant methodological problems, are indeed disrupted in neuropsychiatric disorders and likely have an underlying genetic basis.

ZNF804A Transcriptional Networks in Differentiating Neurons Derived from Induced Pluripotent Stem Cells of Human Origin.

ZNF804A (Zinc Finger Protein 804A) has been identified as a candidate gene for schizophrenia (SZ), autism spectrum disorders (ASD), and bipolar disorder (BD) in replicated genome wide association studies (GWAS) and by copy number variation (CNV) analysis. Although its function has not been well-characterized, ZNF804A contains a C2H2-type zinc-finger domain, suggesting that it has DNA binding properties, and consequently, a role in regulating gene expression. To further explore the role of ZNF804A on gene expression and its downstream targets, we used a gene knockdown (KD) approach to reduce its expression in neural progenitor cells (NPCs) derived from induced pluripotent stem cells (iPSCs). KD was accomplished by RNA interference (RNAi) using lentiviral particles containing shRNAs that target ZNF804A mRNA. Stable transduced NPC lines were generated after puromycin selection. A control cell line expressing a random (scrambled) shRNA was also generated. Neuronal differentiation was induced, RNA was harvested after 14 days and transcriptome analysis was carried out using RNA-seq. 1815 genes were found to be differentially expressed at a nominally significant level (p<0.05); 809 decreased in expression in the KD samples, while 1106 increased. Of these, 370 achieved genome wide significance (FDR<0.05); 125 were lower in the KD samples, 245 were higher. Pathway analysis showed that genes involved in interferon-signaling were enriched among those that were down-regulated in the KD samples. Correspondingly, ZNF804A KD was found to affect interferon-alpha 2 (IFNA2)-mediated gene expression. The findings suggest that ZNF804A may affect a differentiating neuron's response to inflammatory cytokines, which is consistent with models of SZ and ASD that support a role for infectious disease, and/or autoimmunity in a subgroup of patients.

Comparison of REST cistromes across human cell types reveals common and context-specific functions.

Recent studies have shown that the transcriptional functions of REST are much broader than repressing neuronal genes in non-neuronal systems. Whether REST occupies similar chromatin regions in different cell types and how it interacts with other transcriptional regulators to execute its functions in a context-dependent manner has not been adequately investigated. We have applied ChIP-seq analysis to identify the REST cistrome in human CD4+ T cells and compared it with published data from 15 other cell types. We found that REST cistromes were distinct among cell types, with REST binding to several tumor suppressors specifically in cancer cells, whereas 7% of the REST peaks in non-neuronal cells were ubiquitously called and <25% were identified for ≥ 5 cell types. Nevertheless, using a quantitative metric directly comparing raw ChIP-seq signals, we found the majority (∼80%) was shared by ≥ 2 cell types. Integration with RNA-seq data showed that REST binding was generally correlated with low gene expression. Close examination revealed that multiple contexts were correlated with reduced expression of REST targets, e.g., the presence of a cognate RE1 motif and cellular specificity of REST binding. These contexts were shown to play a role in differential corepressor recruitment. Furthermore, transcriptional outcome was highly influenced by REST cofactors, e.g., SIN3 and EZH2 co-occupancy marked higher and lower expression of REST targets, respectively. Unexpectedly, the REST cistrome in differentiated neurons exhibited unique features not observed in non-neuronal cells, e.g., the lack of RE1 motifs and an association with active gene expression. Finally, our analysis demonstrated how REST could differentially regulate a transcription network constituted of miRNAs, REST complex and neuronal factors. Overall, our findings of contexts playing critical roles in REST occupancy and regulatory outcome provide insights into the molecular interactions underlying REST's diverse functions, and point to novel roles of REST in differentiated neurons.

Heat shock alters the expression of schizophrenia and autism candidate genes in an induced pluripotent stem cell model of the human telencephalon.

Schizophrenia (SZ) and autism spectrum disorders (ASD) are highly heritable neuropsychiatric disorders, although environmental factors, such as maternal immune activation (MIA), play a role as well. Cytokines mediate the effects of MIA on neurogenesis and behavior in animal models. However, MIA stimulators can also induce a febrile reaction, which could have independent effects on neurogenesis through heat shock (HS)-regulated cellular stress pathways. However, this has not been well-studied. To help understand the role of fever in MIA, we used a recently described model of human brain development in which induced pluripotent stem cells (iPSCs) differentiate into 3-dimensional neuronal aggregates that resemble a first trimester telencephalon. RNA-seq was carried out on aggregates that were heat shocked at 39°C for 24 hours, along with their control partners maintained at 37°C. 186 genes showed significant differences in expression following HS (p<0.05), including known HS-inducible genes, as expected, as well as those coding for NGFR and a number of SZ and ASD candidates, including SMARCA2, DPP10, ARNT2, AHI1 and ZNF804A. The degree to which the expression of these genes decrease or increase during HS is similar to that found in copy loss and copy gain copy number variants (CNVs), although the effects of HS are likely to be transient. The dramatic effect on the expression of some SZ and ASD genes places HS, and perhaps other cellular stressors, into a common conceptual framework with disease-causing genetic variants. The findings also suggest that some candidate genes that are assumed to have a relatively limited impact on SZ and ASD pathogenesis based on a small number of positive genetic findings, such as SMARCA2 and ARNT2, may in fact have a much more substantial role in these disorders - as targets of common environmental stressors.

Patch-clamp recordings and calcium imaging followed by single-cell PCR reveal the developmental profile of 13 genes in iPSC-derived human neurons.

Molecular genetic studies are typically performed on homogenized biological samples, resulting in contamination from non-neuronal cells. To improve expression profiling of neurons we combined patch recordings with single-cell PCR. Two iPSC lines (healthy subject and 22q11.2 deletion) were differentiated into neurons. Patch electrode recordings were performed on 229 human cells from Day-13 to Day-88, followed by capture and single-cell PCR for 13 genes: ACTB, HPRT, vGLUT1, ßTUBIII, COMT, DISC1, GAD1, PAX6, DTNBP1, ERBB4, FOXP1, FOXP2, and GIRK2. Neurons derived from both iPSC lines expressed ßTUBIII, fired action potentials, and experienced spontaneous depolarizations (UP states) ~2 weeks before vGLUT1, GAD1 and GIRK2 appeared. Multisite calcium imaging revealed that these UP states were not synchronized among hESC-H9-derived neurons. The expression of FOXP1, FOXP2 and vGLUT1 was lost after 50 days in culture, in contrast to other continuously expressed genes. When gene expression was combined with electrophysiology, two subsets of genes were apparent; those irrelevant to spontaneous depolarizations (including vGLUT1, GIRK2, FOXP2 and DISC1) and those associated with spontaneous depolarizations (GAD1 and ERBB4). The results demonstrate that in the earliest stages of neuron development, it is useful to combine genetic analysis with physiological characterizations, on a cell-to-cell basis.

Transcriptome comparison of human neurons generated using induced pluripotent stem cells derived from dental pulp and skin fibroblasts.

Induced pluripotent stem cell (iPSC) technology is providing an opportunity to study neuropsychiatric disorders through the capacity to grow patient-specific neurons in vitro. Skin fibroblasts obtained by biopsy have been the most reliable source of cells for reprogramming. However, using other somatic cells obtained by less invasive means would be ideal, especially in children with autism spectrum disorders (ASD) and other neurodevelopmental conditions. In addition to fibroblasts, iPSCs have been developed from cord blood, lymphocytes, hair keratinocytes, and dental pulp from deciduous teeth. Of these, dental pulp would be a good source for neurodevelopmental disorders in children because obtaining material is non-invasive. We investigated its suitability for disease modeling by carrying out gene expression profiling, using RNA-seq, on differentiated neurons derived from iPSCs made from dental pulp extracted from deciduous teeth (T-iPSCs) and fibroblasts (F-iPSCs). This is the first RNA-seq analysis comparing gene expression profiles in neurons derived from iPSCs made from different somatic cells. For the most part, gene expression profiles were quite similar with only 329 genes showing differential expression at a nominally significant p-value (p<0.05), of which 63 remained significant after correcting for genome-wide analysis (FDR <0.05). The most striking difference was the lower level of expression detected for numerous members of the all four HOX gene families in neurons derived from T-iPSCs. In addition, an increased level of expression was seen for several transcription factors expressed in the developing forebrain (FOXP2, OTX1, and LHX2, for example). Overall, pathway analysis revealed that differentially expressed genes that showed higher levels of expression in neurons derived from T-iPSCs were enriched for genes implicated in schizophrenia (SZ). The findings suggest that neurons derived from T-iPSCs are suitable for disease-modeling neuropsychiatric disorder and may have some advantages over those derived from F-iPSCs.

Allele-biased expression in differentiating human neurons: implications for neuropsychiatric disorders.

Stochastic processes and imprinting, along with genetic factors, lead to monoallelic or allele-biased gene expression. Stochastic monoallelic expression fine-tunes information processing in immune cells and the olfactory system, and imprinting plays an important role in development. Recent studies suggest that both stochastic events and imprinting may be more widespread than previously considered. We are interested in allele-biased gene expression occurring in the brain because parent-of-origin effects suggestive of imprinting appear to play a role in the transmission of schizophrenia (SZ) and autism spectrum disorders (ASD) in some families. In addition, allele-biased expression could help explain monozygotic (MZ) twin discordance and reduced penetrance. The ability to study allele-biased expression in human neurons has been transformed with the advent of induced pluripotent stem cell (iPSC) technology and next generation sequencing. Using transcriptome sequencing (RNA-Seq) we identified 801 genes in differentiating neurons that were expressed in an allele-biased manner. These included a number of putative SZ and ASD candidates, such as A2BP1 (RBFOX1), ERBB4, NLGN4X, NRG1, NRG3, NRXN1, and NLGN1. Overall, there was a modest enrichment for SZ and ASD candidate genes among those that showed evidence for allele-biased expression (chi-square, p = 0.02). In addition to helping explain MZ twin discordance and reduced penetrance, the capacity to group many candidate genes affecting a variety of molecular and cellular pathways under a common regulatory process - allele-biased expression - could have therapeutic implications.

RNA-Seq of human neurons derived from iPS cells reveals candidate long non-coding RNAs involved in neurogenesis and neuropsychiatric disorders.

Genome-wide expression analysis using next generation sequencing (RNA-Seq) provides an opportunity for in-depth molecular profiling of fundamental biological processes, such as cellular differentiation and malignant transformation. Differentiating human neurons derived from induced pluripotent stem cells (iPSCs) provide an ideal system for RNA-Seq since defective neurogenesis caused by abnormalities in transcription factors, DNA methylation, and chromatin modifiers lie at the heart of some neuropsychiatric disorders. As a preliminary step towards applying next generation sequencing using neurons derived from patient-specific iPSCs, we have carried out an RNA-Seq analysis on control human neurons. Dramatic changes in the expression of coding genes, long non-coding RNAs (lncRNAs), pseudogenes, and splice isoforms were seen during the transition from pluripotent stem cells to early differentiating neurons. A number of genes that undergo radical changes in expression during this transition include candidates for schizophrenia (SZ), bipolar disorder (BD) and autism spectrum disorders (ASD) that function as transcription factors and chromatin modifiers, such as POU3F2 and ZNF804A, and genes coding for cell adhesion proteins implicated in these conditions including NRXN1 and NLGN1. In addition, a number of novel lncRNAs were found to undergo dramatic changes in expression, one of which is HOTAIRM1, a regulator of several HOXA genes during myelopoiesis. The increase we observed in differentiating neurons suggests a role in neurogenesis as well. Finally, several lncRNAs that map near SNPs associated with SZ in genome wide association studies also increase during neuronal differentiation, suggesting that these novel transcripts may be abnormally regulated in a subgroup of patients.

Development of patient-specific neurons in schizophrenia using induced pluripotent stem cells.

Induced pluripotent stem cell (iPSC) technology has the potential to transform regenerative medicine. It also offers a powerful tool for establishing in vitro models of disease, in particular, for neuropsychiatric disorders where live human neurons are essentially impossible to procure. Using iPSCs derived from three schizophrenia (SZ) patients, one of whom has 22q11.2del (velocardiofacial syndrome; VCFS), the authors developed a culture system to study SZ on a molecular and cellular level. SZ iPSCs were differentiated into functional, primarily glutamatergic neurons that were able to fire action potentials after ∼8 weeks in culture. Early differentiating neurons expressed a number of transcription factors/chromatin remodeling proteins and synaptic proteins relevant to SZ pathogenesis, including ZNF804A, RELN, CNTNAP2, CTNNA2, SMARCA2, and NRXN1. Although a small number of lines were developed in this preliminary study, the SZ line containing 22q11.2del showed a significant delay in the reduction of endogenous OCT4 and NANOG expression that normally occurs during differentiation. Constitutive expression of OCT4 has been observed in Dgcr8-deficient mouse embryonic stem cells (mESCs); DGCR8 maps to the 22q11.2-deleted region. These findings demonstrate that the method of inducing neural differentiation employed is useful for disease modeling in SZ and that the transition of iPSCs with 22q11.2 deletions towards a differentiated state may be marked by subtle changes in expression of pluripotency-associated genes.